Decay of the Loschmidt Echo for quantum states with sub-Planck scale structures 
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Quantum states extended over a large volume in phase space have oscillations from quantum 
interferences in their Wigner distribution on scales smaller than h [W.H. Zurek, Nature 412, 712 
(2001)]. We investigate the influence of those sub-Planck scale structures on the sensitivity to an 
external perturbation of the state's time evolution. While we do find an accelerated decay of the 
Loschmidt Echo for an extended state in comparison to a localized wavepacket, the acceleration is 
described entirely by the classical Lyapunov exponent and hence cannot originate from quantum 
interference. 
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One common interpretation of the Heisenberg uncer- 
tainty principle is that phase-space structures on scales 
smaller than Ti have no observable consequence. The re- 
cent assertion of Zurek that sub-Planck scale struc- 
tures in the Wigner function enhance the sensitivity of 
a quantum state to an external perturbation, therefore, 
came out as particularly intriguing |^ and even contro- 
versial 1^. His argument can be summarized as follows. 
The overlap (squared amplitude of the scalar product) of 
two quantum states ip and "0' is given by the phase-space 
integral of the product of their Wigner functions, 



miP')\^ = i2TThf / drdpW^W^^. (1) 



For an extended quantum state covering a large volume 
h"^ of 2(i-dimensional phase space, the Wigner func- 
tion exhibits oscillations from quantum interferences 
on a scale corresponding to an action 6S ~ h'^ /A^^"^ <^ h. 
These sub-Planck scale oscillations are brought out of 
phase by a shift 6p, 5x with 5p 5x 5S <^ Ti. The shifted 
state Ip' is then nearly orthogonal to tp since /^,^' ~ 0. 
Zurek concludes that sub-Planck structures substantially 
enhance the sensitivity of a quantum state to an external 
perturbation. 

A measure of_this sensitivity is provided by the 
Loschmidt Echo 



M{t) = |(7/;|exp(ii/t)exp(-ii?ot)|V')P 



(2) 



giving the decaying overlap of two wave functions that 
start out identically and evolve under the action of two 
slightly different Hamiltonians i/o and H ~ Hq + Hi. 
(We set h — 1.) One can interpret M(t) as the fidelity 
with which a quantum state can be reconstructed by 
inverting the dynamics with a perturbed Hamiltonian. 
In the context of environment-induced dephasing, M(t) 
measures the decay of quantum interferences in a system 
with few degrees of freedom interacting with an environ- 
ment (with many more degrees of freedom) Q. In this 
case ijj represents the state of the environment, which in 
general extends over a large volume of phase space. This 



motivated Karkuszewski, Jarzynski, and Zurek Q to in- 
vestigate the dependence of M{t) on short-scale struc- 
tures. 

In this paper we study the same problem as in Ref. , 
but arrive at opposite conclusions. Finer and finer struc- 
tures naturally develop in phase space when an initially 
narrow wavepacket ^pQ evolves in time under the influ- 
ence of a chaotic Hamiltonian Hq As in Ref. [Q, 
we observe numerically a more rapid decay of M(t) for 
ip = exjp{—iHQT)'ipQ as the preparation time T is made 
larger and larger, with a saturation at the Ehrenfest time. 
However, we demonstrate that this enhanced decay is 
described entirely by the classical Lyapunov exponent, 
and hence is insensitive to the quantum interference that 
leads to the sub-Planck scale structures in the Wigner 
function. 

In the case of a narrow initial wavepacket, M(t) 
has been calculated semiclassically by Jalabert and 
Pastawski Before discussing extended states with 
short-scale structures, we recapitulate their calculation. 
The time-evolution of a wavepacket centered at Tq is ap- 
proximated by 

V'(r,i)= /"dro^iff (r,ro;t)7/.o(ro), (3) 

iff (r, ro;t) = Cl''' exp[i5f (r, v,,-t) - i^iis^- (4) 

The semiclassical propagator is a sum over classical tra- 
jectories (labelled s) that connect r and Tq in the time t. 
For each s, the partial propagator is expressed in terms 
of the action integral S^{r, Tq; t) along s, a Maslov index 
fis (which will drop out), and the determinant Cs of the 
monodromy matrix. After a stationary phase approxi- 
mation, one gets 



M{t) 



drViff(r,ro;i)*Xfnr,ro;t)| 



(5) 



Squaring the amplitude in Eq. (|^) leads to a double sum 
over classical paths s, s' and a double integration over 
final coordinates r, r'. Accordingly, M{t) splits into di- 
agonal (s = s', r — r') and nondiagonal (s s' or r r') 
contributions. Since quantum phases entirely drop out 
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of the diagonal contribution, its decay is solely deter- 
mined by the classical quantity Cs oc exp(— At). Here 
A is the Lyapunov exponent of the classical chaotic dy- 
namics, which we assume is the same for H and H^. The 
nondiagonal contribution also leads to an exponential de- 
cay, which however originates from the phase difference 
accumulated when travelling along a classical path with 
two different Hamiltonians |g| . The slope F of this decay 
is the golden rule spreading width of an eigenstate of Hq 
over the eigenbasis of H [§,0. Since M(t) is given by 
the sum of these two exponentials, the Lyapunov decay 
will prevail for F > A. 

The Lyapunov decay sensitively depends on the choice 
of an initial narrow wavepacket ij^n . The faster decay 
of M{t) resulting from the increased complexity of the 
initial state can be quantitatively investigated by consid- 
ering prepared states tp — exp(— «iJoT)'0Oi i-e. narrow 
wavepackets that propagate during a time T with the 
Hamiltonian |12|| , thereby developing finer and finer 
structures in phase space as T increases The sta- 

tionary phase approximation to the fidelity then reads 
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(6) 

with the time-dependent Hamiltonian H.^ = Hq for 
T < T and = i? for r > T. 

We can apply the same analysis as in Ref. to the 
time-dependent Hamiltonian. Only the time interval 
(T, t + T) of length t leads to a phase difference between 
K^^ and K^°, because Hr = Hq ioi t < T. Hence the 
nondiagonal contribution to Mx{t), which is entirely due 
to this phase difference, still decays ex exp(— Ft), inde- 
pendent of the preparation time T. We will see below 
that this is in agreement with a fully quantum mechani- 
cal approach according to which the golden rule decay is 
independent of the complexity of the initial state. 

The preparation does however have an effect on the 
diagonal contribution M^\t) to the fidelity. It decays 
oc exp[— A(t + T)] instead of oc exp(— At), provided t, 
T ^ A^^. This is most easily seen from the expression 

M!,^\t)= f dr5]|X^(r,ro;t + r)|2 

X |iff«(r,ro;t + r)|2, (7) 

by following a path from its endpoint r to an intermedi- 
ate point Ti reached after a time t. The time-evolution 
from r to leads to an exponential decrease oc exp(— At) 
as in Ref. |5j. Due to the classical chaoticity of Hq, the 
subsequent evolution from to Tq in a time T brings in 
an additional prefactor exp(— AT). 

The combination of diagonal and nondiagonal contri- 
butions results in the bi-exponential decay (valid for Ft, 
At, AT > 1) 

Mrit) = A{t) exp(-Ft) + B(t) exp[-A(t + T)], (8) 



with prefactors A and B that depend algebraically on 
time. The Lyapunov decay prevails if F > A and 
t > AT/(F — A), while the golden rule decay dominates 
if either F < A or t < Ar/(F - A). In both regimes the 
decay saturates when Mt has reached its minimal value 
l/I, where / is the total accessible volume of phase space 
in units of h"^. In the Lyapunov regime, this saturation 
occurs at t = t£; — T, where tE = A~^ In / is the Ehren- 
fest time. When the preparation time T ^ Ie, we have 
a complete decay within a time A^^ of the fidelity down 
to its minimal value. 

We now present numerical checks of these analytical 
results for the Hamiltonian 

Ho = (7r/2To)5y + {K/2S)Sl ^ 5{t - nro). (9) 

n 

This kicked top model describes a vector spin of mag- 
nitude S undergoing a free precession around the y-axis 
and being periodically perturbed (period tq) by sudden 
rotations around the z-axis over an angle proportional to 
Sz- The time evolution after n periods is given by the 
n-th power of the Floquet operator 

Fo = eyi^[-i{K/2S)Sl] exp[-i(7r/2)S'y]. (10) 

Depending on the kicking strength if, the classical dy- 
namics is regular, partially chaotic, or fully chaotic. We 
perturb the reversed time evolution by a periodic rota- 
tion of constant angle around the x-axis, slightly delayed 
with respect to the kicks in iJo, 

7/i = (^5,^5(t-nro-e). (11) 

n 

The corresponding Floquet operator is T = 
exp(— i<^S'a:)_Fo- We set tq = 1 for ease of notation. We 
took S = 500 (both H and Hq conserve the spin mag- 
nitude, the corresponding phase space being the sphere 
of radius S) and calculated the averaged decay Mt of 
AfT(t = n) = |(V'|(Tt)"Fo"|V')P taken over 100 initial 
states. 

We choose V'o as a Gaussian wavepacket (coherent 
state) centered on a point (0^ ip) in spherical coordi- 
nates. The state is then prepared as V' = exp(— ^i^oT)■^/'o■ 
We can reach the Lyapunov regime by selecting initial 
wavepackets centered in the chaotic region of the mixed 
phase space for the Hamiltonian (^ with kicking strength 
K — 3.9 1^. Fig. 1 gives a clear confirmation of the pre- 
dicted decay oc exp[— A(t + T)] in the Lyapunov regime. 
The additional decay induced by the preparation time T 
can be quantified via the time tc it takes for Mt to reach 
a given threshold Mc 0. We expect 

tc = -X-^lnMc-T, (12) 

provided > 1// = 1/25 and T < -A^^ Inil^. In the 
inset to Fig. 1 we confirm this formula for Mc = 10^^. 
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As expected, tc saturates at the first kick {tc = 1) when 
T ~ — In Mc < tE- Numerical results qualitatively 
similar to those shown in the inset to Fig. 1 were 
obtained in Ref. and interpreted there as the accel- 
erated decay resulting from sub-Planck scale structures. 
The fact that our numerical data is described so well by 
Eq. ( ]T^ ) points to a classical rather than a quantum ori- 
gin of the decay acceleration. Indeed, Eq. ( [l2|) contains 
only the classical Lyapunov exponent as a system depen- 
dent parameter, so that it cannot be sensitive to any fine 
structure in phase space resulting from quantum inter- 
ference. 
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FIG. 1. Decay of the average fidelity Mt for the kicked 
top with parameters (f) = 1.2 x 10"'^, K = 3.9 and for 
preparation times T = Q (circles), 2 (diamonds), 4 (trian- 
gles), and 6 (squares). In each case, the dashed lines give 
the analytical decay Mt = exp[— A(t -\- T)], in the Lya- 
punov regime with A = 0.42. Inset: threshold time at which 
Mritc) ~ Mc = 10~^. The solid line gives the analytical 
behavior tc = -A"^ \nMc - T. 

We next illustrate the independence of Mrit) on the 
preparation time T in the golden rule regime, i.e. at 
larger kicking strength K when A > F pj. As shown 
in Fig. 2, the decay of Mxit) is the same for the four 
different preparation times T = 0, 5, 10, and 20. We 
estimate the Ehrenfest time as t^; w 7, so that increasing 
T further does not increases the complexity of the initial 
state. 

These numerical data give a clear confirmation of the 
semiclassical result (||). Previous investigations have es- 
tablished the existence of five different regimes for the 
decay of M{t) |,|j9|jl^Jl|l , and since only two of them 
are captured by the semiclassical approach used above, 
we now show that short-scale structures do not affect the 
remaining three. The five regimes correspond to different 
decays: 

(i) Parabohc decay, M{t) = 1 - a^t^, with tr^ = 
(V'ol-ffi iV'o) — (V'ol-ffilV'o)^, which exists for any pertur- 
bation strength at short enough times. 



(ii) Gaussian decay, M{t) oc exp(— cr^t^), valid if a is 
much smaller than the level spacing A. 

(iii) Golden rule decay, M{t) oc exp(— Ft), with F ~ 

ctVa, ff a < f < a. 

(iv) Lyapunov decay, M{t) oc exp(— At), if A < F. 

(v) Gaussian decay, M{t) oc exp(— i?^t^), if Hi is so large 
thxdvat F is larger than the energy bandwidth B of H. 
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FIG. 2. Decay of Mt in the golden rule regime for 
(p = 2.6 X 10~*, 3.8 X 10-"*, 5 x lO""*, K = 13.1, and for 
preparation times T = 0, 5, 10, and 20 (nearly indistinguish- 
able dashed lines). The solid lines give the corresponding 
golden rule decay with F = 0.84 0^5^ ^. 

All these regimes except regime (iii) can be dealt with 
quantum mechanically under the sole assumption that 
both Hq and H are classically chaotic, using Random 
Matrix Theory (RMT) @. Both sets of eigenstates la) 
of H (with N eigenvalues £„) and |ao) of Hq (with N 
eigenvalues e^) are then rotationally invariant [|^. Ex- 
panding ip = X^aV'alo!) and assuming unbroken time- 
reversal symmetry, the fidelity (|^) can be rewritten as 

M{t) = J2 V'aV'/3V'7V'5(a|exp(-zi/ot)|/3) 

X (7I exp{iHot)\S) exp[i{ea - es)t]- (13) 



RMT implies the ^/'-independent average ■(/'a "0/3 ■07 "^iS = 
(SafsS^yS + SofySps + SasSpy) / N"^ . The third contraction 
gives a contribution representing the saturation of 
M{t) for t — > 00. The other two give the time dependence 



M{t) = N-^ + 2N' 



a/3o 



|(a|/?o)pexp[z(e„-e")t]| 



(14) 



For perturbatively weak Hi one has Eq = ej^ -I- (a|-ffi|a) 
and (q;|/3o) — 5a,po- According to RMT the matrix 
elements (a|i?i|a) are independent random numbers, 
and for large N the central limit theorem leads to the 
Gaussian decay (ii) (or the parabolic decay (i) for short 
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times). At larger perturbation strength, |(a|/3o)P be- 
comes Lorentzian, 



(«l/3o)P 



r/27r 



r2/4' 



(15) 



with a width F ~ |(ao|i?i|/3)p/A given by the golden 
rule. This leads to regime (iii). Increasing Hi further 

N-^ and 
This 



one obtains an ergodic distribution |(a|/3o)P ~ 
a straightforward calculation produces regime (v) 
establishes that, under the sole assumption that Hq and 
H are classically chaotic, the decay of the fidelity in the 
three quantum regimes (ii), (iii), and (v) does not depend 
on the choice of the initial state ip. 

In summary, we have investigated the decay of the 
Loschmidt Echo, Eq. (||), for quantum states ^ 
exp{—iHQT)'>pQ that have spread over phase space for a 
time T. As in Ref. we found a faster decay of Mrit) 
than for a localized wavepacket, but only in the regime 
where the decay rate is set by the classical Lyapunov 
exponent A. Since quantum interferences play no role 
in this regime, we conclude that sub-Planck scale struc- 
tures in the Wigner representation of ijj do not influence 
the decay of the Loschmidt Echo. 
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